Lattice stretching bistability and dynamic heterogeneity 
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A simple one-dimensional lattice model is suggested to describe the experimentally observed 
plateau in force-stretching diagrams for some macromolecules. This chain model involves the nearest- 
neighbor interaction of a Morse-like potential (required to have a saturation branch) and an harmonic 
second-neighbor coupling. Under an external stretching applied to the chain ends, the intersite 
Morse-like potential results in the appearance of a double- well potential within each chain monomer, 
whereas the interaction between the second neighbors provides a homogeneous bistable (degenerate) 
ground state, at least within a certain part of the chain. As a result, different conformational changes 
occur in the chain under the external forcing. The transition regions between these conformations are 
described as topological solitons. With a strong second-neighbor interaction, the solitons describe 
the transition between the bistable ground states. However, the key point of the model is the 
appearance of a heterogenous structure, when the second-neighbor coupling is sufficiently weak. In 
this case, a part of the chain has short bonds with a single- well potential, whereas the complementary 
part admits strongly stretched bonds with a double-well potential. This case allows us to explain 
the existence of a plateau in the force-stretching diagram for DNA and alpha-helix protein. Finally, 
the soliton dynamics are studied in detail. 

PACS numbers: 05.45.Yv, 63.20.Ry, 45.90.+t 



> 

m 
oo 
oo 
in 



I. INTRODUCTION 

The effect of appearance of bistable states caused by 
lattice stretching has been first studied by Manevitch et 
al. [l[ for modeling the mechanodestruction of a polymer 
chain. As a simple model, the anharmonic chain with the 
nearest-neighbor coupling in the form of the Morse po- 
tential has been chosen. Under lengthening this chain by 
applying an external force to its ends, the formation of 
an effective double-well potential in the chain bonds has 
been shown. The idea of this conformational transition 
can be explained in simple terms as follows. Consider 
three coupled particles as shown in Fig. [TJ where the two 
lateral particles are fixed and the central particle inter- 
acts with its neighbors through a Morse-like potential 
V(r) with a minimum at R = tq and a constant asymp- 
tote e = linv^oo V(r). The potential of this type has a 
point of inflection at r = Rq > tq. The Morse potential 



form as follows 



V(r) = e 



-0(r-r o ) 



< r < oo, 



(1) 



given in dimensionless units, with any parameter /3 > 0, 
can be chosen as a particular example where Rq = tq + 
/3 _1 ln2. The total potential for the middle particle is 
V(r) + V(2R — r), where 2R is the distance between the 
lateral particles. Using the new variable u = R — r, this 
potential can be written in a more convenient symmetric 



-2f3{R-r ) 



W{u) = V(R -u) + V(R + u) = ee 
x{[2 cosh(/3u) - e ^ R - ro) ] 2 + e 20{R - ro) - 2}. (2) 

The potential W(u) has only one minimum u = if 
R < Rq as demonstrated by Fig. DJa), and one max- 
imum u = and the two minima u — ±uo if R > 
R [see Fig. QJb-d)]. For the potential Q, u {R) = 
/3~ 1 arccosh(e' 3 ( fl-ro ' ) /2) > 0, being the solution of the 
equation 



cosh(/3u ) 



_ „/3(fl-r ) 



/2- 



(3) 



The three-particle system illustrated by Fig. Q] can be 
extended to a finite chain consisting of 2N + 1 particles 
(or N monomers), where the terminal particles are fixed 
and the total chain length is L = 2NR. If R > Rq, 
double- well potentials can be formed inside the chain. In 
this case, many ground states of the chain are possible re- 
sulting in different irregular chain conformations. There- 
fore the model studied in Ref. [l[ has to be modified in 
such a way that a sufficiently stretched chain would ad- 
mit homogeneous ground states with periodic structure. 
To this end, we involve additionally a stabilizing second- 
neighbor interaction and, as a result, the ground states 
with alternating lengths of chain bonds are possible. 

It is sufficient to impose an harmonic coupling be- 
tween the second neighbors as shown schematically in 
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FIG. 1: (Color online) (a) Single-well and (b-d) double- 
well potentials for a central particle interacting with its 
fixed neighbors via inter-atomic Morse-like potential V(r) and 
formed as the sum V(r) + V(2R— r). (a) Distance 2R between 
the lateral fixed particles does not exceed 2Ro, where Ro is a 
point of inflection of V(r). (b) Unstable equilibrium position 
of the central particle, (c) Left and (d) right stable equilibria 
of the central particle. 



Fig. [21 Let K be a (dimensionless) stiffness constant of 
the second-neighbor interaction with xq, x\, . . . , x-2N be- 
ing positions of the chain atoms. Then the total poten- 
tial energy of the TV-monomer system with fixed terminal 
atoms can be written as 



2N-1 



En — 



~ ^2 ^(zn+i - x n ) + K(xi - x - r Q ) 2 

71 = 

2N-2 ^ 

+ ^ n K ( X n+2 - X n - 2r ) 2 



n=0 



-K(x 2 n - X2N-X - r ) 2 



(4) 



It is expected that a sufficiently strong stretching of the 
chain results in a dimcralization of the chain, for which 
the even atoms are found in equilibria x n — nR with 
n = 0,2,..., 2N. Inserting these values into the energy 
Q, we find 



E N = N [V(r) + V(2R - r)} + 4(/V - 1)K(R ~ r ) 2 

r-r ) 2 . (5) 



+ \K{r - r ) 2 + \K{2R- 



Differentiating Eq. ([5]) with respect to r, we find the fol- 
lowing equation for equilibria: 



V'(r) - V'(2R - r) = {2K/N)(R - r). 



(0) 



Using the variable u = R — r, Eq. (|6|) becomes W'(u) — 
—2Ku/N. The trivial solution u = describes equilibria 
of odd atoms (stable if R < Rq and unstable if R > Ro). 
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FIG. 2: (Color online) Equilibria in monoatomic chain with 
fixed terminal particles, where except for the nearest-neighbor 
interactions, also a coupling between the second neighbors is 
involved: (a) Stable (ground) state of the stretched chain with 
period R < Ro, where each atom is found in a single- well po- 
tential, (b) Unstable state of the chain, where all atoms with 
odd number are located at the top of a double-well poten- 
tial, (c) Left and (d) right degenerate ground states of the 
stretched chain. 



The two stable equilibria with u — ±uq appear when 
R > Rq and N » 1. For a long chain (N — > oo), the role 
of boundary conditions can be neglected, so that Eq. (|6]) 
for the equilibrium positions ±ito in the case of the Morse 
potential ([1]) takes the form of Eq. ([3]), so that the two 
stable minima exist if the inequality R > Ro is fulfilled. 

To simplify the problem with fixed chain ends, it is con- 
venient to use the cyclic boundary conditions by putting 



X2N 



x + 2NR and x 2N +\ =xi + 2NR. (7) 



Then, if we fix the whole chain by setting xo = 0, so that 
x\ = ±mq, the equilibrium positions are 



nR if n 

nR ± uq if n - 



0,2,. 
1,3,.. 



,2iV, 
,27V- 



1, 



(8) 



The subscripts " + " and "-" denote the two degen- 
erate ground states in the dimeralized chain, respec- 
tively. Schematically, these states can be represented as 

X-X X-X X ■■■ X-X X-X X and 

X X-X X-X ■■■ X X-X X-X, where 

the terminal X's are fixed and all the bulk atoms are 
found either in the left well or the right well, respec- 
tively. Obviously, the domain walls (topological kinks 
and antikinks) that separate these two ground states can 
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be excited. However, this is true if the second-neighbor 
interaction is sufficiently strong, at least for the model 
suggested in this paper. As shown below, the situation 
appears more complicated for a weak second-neighbor 
coupling, the case being of experimental relevance for 
some macromolecules. More precisely, the existence of a 
plateau in the force-stretching diagrams for DNA double 
helix [2 6] as well as for a-helices of protein [7[ can be 
explained within the framework of our model. 

The paper is organized as follows. In the next sec- 
tion, we present the equations of motion for a stretched 
nonlinear monoatomic chain and discuss the spectrum of 
small-amplitude oscillations. The analysis of switching 
a bistable ground state of the chain is given in Sec. III. 
In Sec. IV, we find the profiles of kink and antikink so- 
lutions and study their dynamical properties. The next 
section is devoted to realistic systems, where the topo- 
logical soliton solutions obtained in the previous section 
are studied in detail. Conclusions are given in Sec. VI. 



II. A MODEL AND ITS LINEARIZED VERSION 

With the notations introduced in the previous section, 
the (dimensionless) Hamiltonian for the monoatomic 
chain model with the cyclic boundary conditions ([7} can 
be written in the form 



2JV-1 r 



^mxl + V(x n+1 - x n ) 



+ -K(x n+2 - x n - 2r ) 2 



where m is a chain particle mass and the dot denotes the 
differentiation over time t. Here the strings connecting 
the second neighbors are assumed to be undistorted at 
the length 2r$. The corresponding equations of motion 



mx n = [V'(x n+ i -x n )- V'{x n - x n -i)\ 

+K(x n+2 -2x n + x n ^ 2 ), (10) 

where n = 0, 1, . . . , 2iV - 1. 

The linearized version of Eqs. (fl"0| is obtained by 
putting x n (t) = £° + v n (t) , n = 0, 1, ... , 2N + 1, where 
the equilibria defined by Eqs. (jU). As a result, 

we find 



[k±(u 



n+l 



+K(v n+2 - 2v n + v n - 2 ), 



(11) 
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FIG. 3: (Color online) Stiffness k — V" (R±uo) as a function 
of lattice constant R calculated for potential Q with e = 1, 
/3 = 1, and ro = 1: k = k± for R < Ro (curve 1), k — fc+ 
(curve 2) and k — k~ (curve 3) for R > Rq. Critical stiffness 
k — k c = — fc_fc+/2(fc_ + k+) is given by curve 4. 



k = k- = k + = V" (R) . With the stretching of the chain, 
the stiffness k(R) decreases monotonically and at R = Ro 
it reaches zero (Fig. |3l curve I). At further lengthening, 
the chain becomes bistable and the distance 2uq increases 
monotonically with the growth of R. The stiffness of the 
short bonds k- = eV"(R — uo) > increases monotoni- 
cally (Fig. |31 curve 3) and the stiffness of the long bonds 

(9) 

becomes negative: k+ — V"(R + u ) < (curve 2), but 



their sum is always positive: fe- 
ttle potential ((T|) we get 



-k + > 0. Explicitly, for 



k± = Te/3 z e TPu ° sinh(/?u ) costT z {/3u ), 



(12) 



so that in this case fe_ + k + = 2ef? tanh 2 (/?u ) > 0. 

The dispersion law is obtained if we insert the small- 
amplitude waves 

v n (t) = Ae i(qnR - ut} and v n+1 (t) = Be i ( qnR ~ ut) (13) 

into Eqs. flTTJ with n = 0, 2, . . . , 2N . Here u is the fre- 
quency and the wave number q € [0, n/2R). As a result, 
we find the dispersion law that admits two branches of 
the spectrum (acoustic and optical): 



fc_ + fc+ + 2K[1 - cos(2qR)} 



±Jk 2 _+kl + 2k-k+cos(2qR). (14) 



where k± = V"(R±uq). The upper subscript at the stiff- 
ness coefficient belongs to the particles with even (odd) 
n's and the lower one to those with odd (even) n's for 
the case when all the odd particles are found in the right 
(left) well of the double-well potential. 

Under stretching R < R , the displacement u be- 
comes zero and therefore we have the stiffness constant 



The dispersion curves for R < Rq are present in 
Fig. HJa). For R < Rq the curve 0J-(q) is the con- 
tinuation of the curve ui+(q). They can be considered 
as a single acoustic branch. In the limit R — » Rq — 0, 
the curves ui-(q) and w+(g) merge together into a single 
curve uJ 2 = 2K[1 — cos(2gi?)]. At further increase of R, 
this curve splits into two disconnected curves w_ (q) and 
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2qR 



where the dependence of u$ on R is given by 
Eq. ([3]). Next, we find the limiting value Kq = 
UmR-^.R 0+ o k c (R) = e/3 2 /4 and therefore the stability 
condition of the alternating states of the stretched chain 
takes the following simple form: 

K > K . (18) 

Thus, for the stability of the stretched alternating states 
of the chain, it is necessary and sufficient that the stiff- 
ness of the interaction of the second neighbors has to be 
greater the eighth part of the stiffness of the interaction 
the nearest neighbors of the unstretched chain. 

The velocity of long-wave acoustic phonons vq = 
lim 9 ^ u{q)/q can be calculated from the spectra (fT4")l 
for different values of R. As a result, we obtain 

_ / 2fly W4 + g)/m if R<Rp, ( 
\ 2Ry/(K-k c )/m if R > R Q . 1 ; 

Note that the condition (fTS")) ensures the positivity of the 
expression under the radicals in ([T9|) . 



FIG. 4: (Color online) Linear spectra cj — oj±(2qR) calculated 
(at m = 1) in the case of potential {TJ with e = 1, /? = 1, 
ro = 1, and K = 0.5 for (a) R = ro (curves 1 and 2), R = Ro 
(curve 3), and (b) R = 2 > Ro (curves 4 and 5). Curves 
1,3, and 4 relate to uj-(2qR) and curves 2,3, and 5 relate to 
UJ+{2qR). 

uj+(q) [Fig. QJb), curves 4 and 5], where the first curve 
corresponds to acoustic oscillations, while the second one 
to optical oscillations of the stretched chain. 

At stretching R < R®, the uniformly stretched state 
is always stable, since all the coupling constants fc_, fc+, 
and K are positive. The situation changes when R > 
Ro because here one of the coupling constants, fc+, is 
negative and now the stability of the alternating state of 
the stretched chain depends on the stiffness constant K . 
For the stability it is necessary that the inequality 

k- + k+ + 2K[1 - cos(2qR)} 

-\Jk 2 _+k\ + 2k^k+ cos(2qR) > (15) 

has to be fulfilled for all values of the wave number q > 0. 
It is easy to show that this condition holds only if 

K > k c {R) = -k_k + /2(k_ + k + ). (16) 

The dependence of the critical value of the stiffness of 
the second-neighbor interaction on R is given in Fig. [3] 
(curve 4). For all R this value is positive; it monoton- 
ically decreases with increase of the chain stretching R. 
Its maximum reaches in the limit R — > Rq + when 
u —> u + 0. Using Eq. (fT2)) . for the particular case (fTJ) 
we obtain 

k c (R) = £^! C osh- 2 (/3u ), (17) 



III. TRANSITION TO THE BISTABILITY OF 
THE GROUND STATE UNDER STRETCHING 
THE CHAIN 

In order to understand how the ground state of the 
chain changes under its stretching, we consider the de- 
pendence of the ground energy E of the homogeneous 
chain state on the lattice spacing R. For the uniformly 
stretched chain state, when x n +i — x n — R, < R < oo, 
n = 0,1,..., 2N, the deformation energy of one chain 
unit is 

Ej.{R) = V(R) + 2K(R~r ) 2 

= e[e- 0{R - ro) - l} 2 + 2K(R-r ) 2 . (20) 

On the other hand, when we consider the ground state 
with the alternating bond lengths R—uq and R+uq (i? > 
i?o)j the deformation energy of one chain unit becomes 

E 2 (R) = W(u )/2 + 2K(R-r ) 2 

= ^[l-2e' 2 ^ R - ro) ]+2K{R-r ) 2 . (21) 

For comparison the form of the functions Ei(R) and 
E 2 {R) is depicted in Fig. [5] The function Ei(R) has a 
minimum at R = ro, increasing for R > r$. At R = Rq 
both these functions are smoothly "sewed" together be- 
cause £i(i? ) = E 2 (R ) = e/4 + 2K(R - r ) 2 and 
E[(Ro) = E' 2 {R Q ) = e/3/2 + 4K{R - r ). However, 
for R > Ro the function E 2 {R) steps aside smoothly 
and continues further below Ei(R). Therefore the en- 
ergy of the homogeneously stretched (with any R) ground 
state of the chain is given by the smoothly sewed func- 
tion E(R) = E^R) for R < R and E(R) = E 2 (R) for 
R > Rq. At R = R the second derivative E'^Ro) = 
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AK — ef3 2 is positive if K > Kq and negative if K < K$. 
Therefore E(R) is a strongly concave function (for all R) 
only if K > Kq [see Fig. H5T b)]. In this case, the ground 
state of the chain is always the homogeneous conforma- 
tion with equal bond lengths for R < Rq and that with 
alternating bonds for R > Rq. The inequality (fT8j) en- 
sures the stability of the uniformly stretched state of the 
chain. 

For K < Kq a local convexity in the E(R) behav- 
ior, as illustrated in Fig. [IJa) by curve 2, appears in a 
neighborhood of Rq, i.e., on some interval Ri < R < R2 
with Ri < Rq and R2 > Ro- This means that the 
homogeneous state given by the energy E 2 (R), R\ < 
R < R 2 , with the alternating bond lengths R — uq(R) 
and R + uq(R) in fact is unstable. Instead, a heteroge- 
neous conformation, where some part of the chain has 
equal bonds and the other one alternating bonds, ap- 
pears more stable. In this case, the chain energy be- 
havior can be obtained by connecting the two points 
{Ri,E(Rx)} and {R 2 ,E(R 2 )} by a line [see Fig. G^a), 
line 3]. In other words, moving along this line, the het- 
erogeneous state with one part of the chain being in a 
weakly stretched state with equal bonds and the spacing 
i?i , and the other part in a strongly stretched state with 
alternating bond lengths R 2 — uq(R 2 ), R 2 + uq(R 2 ) and 
the spacing R 2 appears more energetically favorable. In 
this case, the heterogeneous stretching (lengthening) of 
the whole chain occurs due to the increase of the por- 
tion of strongly stretched bonds. This scenario of the 
heterogeneous stretching results in the appearance of the 
stationary region (plateau) in the force-stretching dia- 
gram under the chain lengthening < R < R 2 ) as 
illustrated by line 6 in Fig.(5jc). 



IV. DYNAMIC HETEROGENEITY AND 
TOPOLOGICAL SOLITONS 
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FIG. 5: (Color online) Dependence of energy E and its deriva- 
tive dE/dR of uniformly stretched state of the chain with 
equal bonds E = E\ (curves 1 and 4) and that with alternat- 
ing bonds E — E2 (curves 2 and 5) on lattice spacing R for 
(a) and (c), line 6, K = 0.05; and (b) and (c), line 7, K = 0.3 
(/3 = 1, e = 1, r = 1, Ko = 0.25, and R = l + ln2 = 1.693). 
Lines 3 and 6 give convex neighborhood of function E(R). 



Since the nonlinear lattice model introduced in the pre- 
vious section admits the heterogeneous structure that ap- 
pears to be energetically favorable, the existence of freely 
moving topological defects is expected. The correspond- 
ing soliton solutions can be found numerically using the 
steepest-descent method. To this end, it is convenient 
to use the variables: coordinates u n — x n /ro, energy 
H = H/e, and time r = r^ 1 yje/mt. Then the dimen- 
sionless Hamiltonian of the cyclic chain takes the form 



2JV-1 



71=0 



1 / d,U n 
2\dT 

1 



+ V(u n+ i - U n ) 



+ ^K(u n+2 - U n - 2) 2 



(22) 



where V(tt n +i — u n ) = e^ 1 V{x n+ i —x n ) and k — Kr\je. 
According to Eqs. (O, the chain tension is given through 
the boundary conditions: 

u 2 n = uq + 2Na and = = u\ + 2iVa, (23) 



where the lattice spacing of the stretched chain is a — 
R/r > 1. 

The system of equations 



dr 2 



V'(tt n+ i - U n ) - V(u n - U n -i) 
+ n(u n+2 - 2u n + it n _ 2 ), 



(24) 



where n — 0,1,..., 2N — 1, corresponds to the Hamil- 
tonian function (|2"2"|) with the boundary conditions (|2"3")l 
and U-i = u 2 n-i — 2Na and u- 2 = M2JV-2 — 2Na. For 
the relative displacements r n = u n +i — u n , the equations 
of motion (|2^|) become 



d_Tn 
dr 2 



= V'(r„+i) - 2V'(r n ) + V(r«_i)) 
+ K(r n+2 - 2r n +r„_ 2 ). 



(25) 



For numerics we choose the following values of the 
parameters: e = 1, /3 = 1, and tq = 1. Then for 
k > kq = 1/4, the ground state will always be uni- 
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formly stretched because the stability condition of the al- 
ternating states (fP5)l is fulfilled. Consider the chain with 
the lattice spacing a > a — Rq/tq — 1 + (/3r ) _1 ln2 
and k > kq. Then the chain has the following two 
ground states with equal energy: r 2n _ \ — a =p 5 and 
r 2 „ = a ± 6, n — 1, . . . ,N/2, where S — uq/tq. We 
look for traveling wave solutions of the equations of mo- 
tion (|25p that describe the smooth transition of the chain 
from one ground state to the other one, i.e., we put 

rin-iir) = r 1 (2na - sr), r 2n (r) = r 2 (2na - sr), (26) 

where s is a dimensionless traveling wave velocity. Sup- 
pose next that the form of the traveling wave ([26]) 
smoothly depends on the lattice number n. Then the 
second derivatives over time can approximately be sub- 
stituted by the discrete derivatives as follows: 



dr 



dr 2 



r-2,, 



s 2 {r 2n -3 - 2r 2n -i + r 2 „+i)/4a 2 
s 2 {r 2n -2 - 2r 2n + r 2n+2 )/4a 2 , 



(27) 



so that the equations of motion (|25[) transform to the 
system of discrete equations 



V(r»+i) - 2V'(r n ) + V(r„_i) 
+ k(1 - s 2 )(r„ +2 - 2r n + r„_ 2 ) = 0, 



(28) 



where s — s/2a^/n is the reduced value of the velocity. 
It is convenient to look for the solution of the system of 
discrete equations numerically as a solution of the condi- 
tional minimum problem: 

N 1 
F = £){V(r2„-i) + V(r 2n ) + -k(1 - s 2 ) 

n=l 

x [(ran-i + r 2n - 2) 2 + (r 2 „ + r 2n+1 - 2) 2 ]} 

N 



min : ^(r 2 „_ 



+ ra„) = 2Na. (29) 



The problem l|29p has been solved numerically by the 
method of conjugated gradients described in Ref. [8(. The 
initial point corresponding to the presence of the kink- 
antikink pair with the centers at the sites NJ 4 and 3N/ 4 
has been used as follows: 

r 2n -i = 1-6, 

r 2n = 1 + 5 for n < N/4, n > 37V/4, 

r 2n -i = 1 + 5, 

r 2n = 1-8 for N/4 < n < 3N/4. (30) 

Let {r n } 2 l N 1 be a soliton solution of the problem (|2U1) . 
Then it is possible to find the soliton energy 



1 N 

£=-£{V(r 2 „_x) 



V(r 2r , 



x [(ran-i + ^2„ - 2) 2 + (ra„ + r 2n+1 - 2) 2 ]}, (31) 
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FIG. 6: (Color online) Stationary topological soliton (s = 0) 
in stretched chain at a = 1.7, /3 = 1, ro = 1, k = 1, N = 500. 
Dependence of (a) relative displacements ran-i ?~2n (solid and 
dashed lines) and (b) local compression p n — r2n— l + J*2n on 
the number of the chain site are shown. 



and its diameter 



D = l + 2 



N/2 



n)p n 



1/2 



(32) 



where the soliton center is given by 



N/2 



1 



n=l 



(33) 



and the sequence 



(r 2 



n+l 



■r 2n -i)/S, S 



N/2 

^{r2n+l 



■r 2 „_i), (34) 



determines the distribution of deformation along the 
chain. 

The shape of the topological soliton is presented in 
Fig. [5] The panel (a) shows that the lengths of odd 
bonds r 2n _i have the kink shape, whereas the lengths 
of even bonds r 2 „ the antikink shape (and vice versa for 
the antikink). Next, as shown in the panel (b), the lo- 
cal compression of the chain takes place in the region of 
soliton localization. 

The energy of formation of the kink-antikink pair can 
be defined as the difference AE = E — Eo, where E is the 
energy of the stationary kink-antikink pair in the cyclic 
chain and Eq is the energy of the ground state of the chain 
at a given lattice spacing a > clq- As shown in Fig.JTJa), 
the formation energy monotonically increases with the 
growth of the lattice spacing. Nearby the critical value of 
the lattice spacing clq, the energy of formation becomes 
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FIG. 7: (Color online) Dependence of the energy of formation 
of soliton-antisoliton pair E (a) and the diameter of topologi- 
cal soliton D (b) on the value of lattice spacing a of stretched 
cyclic chain consisting of N = 500 sites (/3 = 1, Vq = 1, k = 1, 
and oo = 1.693). 



FIG. 8: Region of existence of topological solitons in the space 
of parameters s, k (a) and region of parameter values at which 
soliton solutions are absent (b). The line which splits these 
regions is given by (|35[) . Lattice spacing is a — 1.7 and j3 = 1, 
r = 1. 



infinitesimal. When a — s- clq + 0, the energy AE 
and the soliton diameter D — > oo. With increasing a, 
the soliton diameter monotonically decreases down to the 
value D = 3.4 [see Fig.^b)]. 

Consider the dependence of the energy and the diam- 
eter of the soliton on its velocity. We choose the values: 
P = 1, r = 1, a = 1.7 > a = 1 + ln2 = 1.693 and 
k > kq ~ 1/4. It follows from Eq. (IT9l that the reduced 
dimensionless velocity of sound is 



so = s /2ayfn = . 1 



K-K+ 



2k(k- + K + ) ' 



(35) 



where k± — V"(a ± S). The numerical solution of the 
problem ([29)) has shown that the system of discrete equa- 
tions (|28[) has a soliton solution only in the subsonic 
region: s < sq, k > kq, being a typical situation for 
topological solitons (kinks) . For our model, the region of 
the existence of topological solitons in the space of the 
parameters k, s is shown in Fig. [5] The region of the 
existence of solitons is separated from the region of their 
absence by curve (|35|) which determines the dependence 
of the reduced velocity of sound so on the dimensionless 
stiffness k. 

The dependence of the energy and the diameter of the 
topological soliton on its velocity is given in Fig. O As 
follows from this figure, the soliton energy monotonically 
increases with the growth the velocity. The energy tends 
to infinity at the velocity of long-wave acoustic phonons. 
The soliton diameter non-monotonically depends on its 
velocity. For small values the velocity increase results 
in negligible decrease of the diameter, which nearby the 




(a) 



(b) 



FIG. 9: (Color online) Dependence of energy E and diameter 
D on reduced velocity s at ft = 1, ro = 1, n = 1, and a = 1.7. 



right edge of the velocity spectrum turns into the fast 
monotonic growth. 

Consider the dynamics of a kink-antikink pair in a 
cyclic chain consisting of N sites. To this end, we in- 
tegrate the system of the equations of motion (|25|) with 
the initial conditions 

r n (0)=r°, for n = l,2,...,2iV; 
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FIG. 10: Uniform motion of kink-antikink pair in stretched 
chain. Dependence of chain distribution of odd t2n-i (a) 
and 6V6I1 T"2ri (b) bond lengths on time r. Parameter values: 
f3 = 1, ro = 1, k = 1, a = 1.7 (lattice spacing of stretched 
chain), s = s/2ay r K = 0.5 (reduced soliton velocity). 



rL-i(0) = s(r 2n+ i ~ r° 2n _ 3 )/4a, (36) 

for n = 1,2,..., N, 

where s is a soliton velocity and {v°}^=i is a soliton 
solution of the conditional minimum problem (|29[) . 

The numerical integration of the system with the 
initial conditions (|36p has shown that the topological soli- 
tons in the stretched chain are dynamically stable for all 
admissible velocities s < s . As illustrated by Fig. [TU1 
the solitons move along the chain with a constant veloc- 
ity without phonon radiation, completely retaining their 
initial shape. 

Consider now the interaction of the solitons with op- 
posite polarity under their collision. To this end, we in- 



(a) 




FIG. 11: Annihilation of two topological solitons with oppo- 
site polarity under collision in stretched chain. Dependence 
of the chain distribution odd rm-x (a) and even ri n (b) bond 
lengths is shown. Parameter values: ft = 1, ro = 1, K = 1, 
a = 1.7 (lattice spacing of stretched chain), s — s/2ay / re = 0.5 
(reduced soliton velocity). 

tegrate the system with the initial conditions 

r„(0) = r°, for n = 1, 2, . . . , 2AT; 
r'2n-i(0) = -s(r 2n+1 -r 2n _ 3 )/ia, 
r'2n(0) = -s(r 2n+2 ~r° 2n _ 2 )/4a, 

for n = 1,2,..., N/2; (37) 
r' 2n -i(0) = s(r° n+1 -r° 2n _ 3 )/Aa, 
r2n(0) = s(r 2n+2 -r° 2n _ 2 )/4a, 
for n = N/2 + 1, N/2 + 2,...,N. 

The numerical integration has shown that this interaction 
is inelastic. Thus, at the velocity s = 0.5, the collision 
results in the annihilation of solitons with opposite po- 
larity. The collision is accompanied by intensive phonon 
radiation and leads to the appearance a breather-like lo- 
calized oscillation (see Fig. ITT1) . 



FIG. 12: Non-uniform stationary state of stretched cyclic 
chain with a > ao and k < Kq (/3 = 1, ro = 1, a = 1-7, 
k = 0.2). Insets show that the edges of strongly stretched 
region of the chain have soliton shape (bond lengths r2 n -i, 
T2n smoothly depend on n). 



For k < kq the chain can be found in the ground state 
of three types, depending on the lattice spacing a. Un- 
der weak stretching a < a\ (1 < a\ < ao), the ground 
state is a uniformly stretched chain with equal bonds. At 
middle stretching a\ < a < a 2 (a 2 > io): a part of the 
chain is found in a weakly stretched homogeneous state 
with equal bonds and the spacing a = <zi, whereas the 
complementary part in a strongly stretched state with 
alternating bonds and the spacing a — a 2 . For a > a 2 
the whole chain is found in a homogeneous state with 
alternating bonds. 

The numerical analysis confirms the existence of these 
three types of chain states. Thus, at f3 = 1, ro = 1, 
k = 0.2, the critical values of the lattice spacing are ao = 
1 + In 2 = 1.6931, ai = 1.6833, a 2 = 1.8655. At a = 1.7 
the main part of the chain is found in the uniform state 
with equal bond lengths ai < ao, whereas the other part 
turns into the strongly stretched alternating state with 
period a 2 > ao (see Fig. IT2"j) . This behavior of the chain 
under stretching can be explained by the non-convexity 
of the function E(a) (see Sec. III). 

Figure [T2l also illustrates that the edges of the strongly 
stretched region of the chain with the alternating struc- 
ture have the form of smooth stairs describing a smooth 
transition of the chain from the state with equal bonds to 
the state with alternating weakly and strongly stretched 
bonds. The numerical simulations have shown that the 
strongly stretched region can propagate along the chain 
with a subsonic velocity completely retaining its shape 
(see Fig.[T3|). Therefore this transition region of the chain 
is a soliton describing the transition of only a part of 
the chain into the state with alternating bonds. As fol- 
lows from this Figure, the collision of these solitons does 
not result in their destruction, but it is accompanied by 
phonon radiation. 
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FIG. 13: Inelastic collision of two strongly stretched regions 
in cyclic chain with a > ao and « < kq (/9 = 1, ro = 1, 
a — 1.7, k — 0.2, velocity is s = 0.25). Dependence of chain 
distribution of odd T2n-i (a) and even vm (b) bond lengths 
on time r is present. 



V. EXAMPLES OF MOLECULAR CHAINS 
WHERE THE EXISTENCE OF STRETCHING 
SOLITONS IS POSSIBLE 

We have studied the simplified one-dimensional lattice 
model. Nevertheless, this study allows us to define a 
family of molecular systems with quasi-one-dimensional 
structure in which the existence of topological solitons of 
stretching is possible. 

The stability condition of alternating states of the 
stretched chain (|18p imposes an important constraint. If 
we suppose that the stiffness of interaction of molecules 
is proportional to the energy of their interaction, then 
this condition can be treated as the requirement that 
the energy of interaction of the second neighbors in a 
quasi-one-dimensional molecular chain has to be no less 
than eighth part of the energy of interaction of the near- 
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est neighbors. This condition of commensurability of the 
energies of interaction cannot obviously be realized in the 
chains where the nearest neighbors are coupled by strong 
valent bonds while the second neighbors by weak Van der 
Waals interactions. 

Note that the condition (TT5|) provides the stability of 
a stretched chain only with respect to its longitudinal 
deformations. In the three-dimensional space, the insta- 
bility of a stretched chain can be caused by other (orien- 
tational, bending, twisting, etc.) deformations. There- 
fore the commensurability condition in this case is not 
sufficient. 

Consider a zigzag-like chain of hydrogen bonds as the 
first example. Hydrogen fluorides HF, chlorides HC1, 
HBr, and HI (HX) at low temperatures have a crystalline 
structure formed by planar zigzag-like chains of hydro- 
gen bonds Consider an isolated hydrogen-bonded 
chain (HX- ■ - )oo consisting of two-atom molecules of flu- 
oride HF and chloride HC1. 

The interaction of two-atom polar molecules HX is usu- 
ally described by the 12-6-1 potential [l3[ 

tl = l 12=1 

with the seven free parameters: two Lennard- Jones pa- 
rameters e and a; three charges q\, q 2 , and q 3 (qi + q 2 + 
<73 = 0), lying on the line of valent bonds, and three 
distances r%, r 2 , and r 3 which assign the charge posi- 
tions. Here is the distance between the charge 
of the first molecule HX and the charge qi 2 of the second 
molecule given in terms of r l5 r 2 , and r 3 . 

The values of the parameters for potential (I3"5)) can 
easily be found using the data of the crystalline structure 
of (HX)^ and ab initio calculations of the dimer (HX) 2 
[l|. We get 

qi = -0.6397e, q 2 = 0.6159e, q 3 = 0.0238e, 
n = 0.25 A, r 2 = 0.9075 A, r 3 = -1.575 A, 
e = 0.00798 eV, a = 2.837 A (39) 

for hydrogen fluoride HF and 

qi = -0.3147e, q 2 = -0.7974e, q 3 = 0.4827e, 
ri = 1.296 A, r 2 = -0.275 A, r 3 = -0.838 A, 
e = 0.0298 eV, a = 3.602 A (40) 

for hydrogen chloride HC1, where e is the electron charge. 

In the plane of the zigzag-like chain (■ • -HX- • ^oo, the 
position of each molecule HX is given by the coordi- 
nates x and y of the center of the heavy molecule X 
and the angle 4> which shows the direction (orienta- 
tion) of the molecule HX. The detailed description of 
the quantum-mechanical model of this zigzag-like struc- 
ture is given in (l5j . For hydrogen fluoride, the zigzag 
angle is a = 119.5°, the longitudinal lattice spacing 
l x = 2.167 A, the distance between the nearest molecules 



po = 2.509 A, the direction of each chain molecule dif- 
fers from the zigzag line only by the angle <po = 01.21°. 
(The parameters for hydrogen chloride are a — 93.6°, 
l x = 2.692 A, po = 3.694 A, tp = 0.95°). In equilib- 
rium, the energy of interaction of the nearest molecules 
is Ei = 0.2339 eV and the energy of interaction of 
the second-neighboring molecules is E 2 = 0.0312 eV 
(£1 = 0.0890 eV, E 2 = 0.0165 eV). 

Here the commensurability condition (p~8j) of the en- 
ergy of interaction of the first and second neighbors is 
fulfilled. Thus, for the chain of molecules of hydrogen 
fluoride E 2 jE\ = 0.133 > 1/8 and for hydrogen chloride 
E 2 jE\ = 0.185 > 1/8. However, the analysis of spec- 
trum behavior under the chain stretching has shown that 
the stability of a uniformly stretched state of the chain 
disappears before reaching the point of inflection of the 
effective potential of longitudinal stretching. In the case 
of HF, for the point of inflection the longitudinal lattice 
spacing of the zigzag is ao = 2.67 A, while the stability of 
the chain disappears already at the longitudinal lattice 
spacing a x = 2.625 A (for HC1 a = 3.60 A, a x = 3.48 A). 
For the longitudinal spacing of the zigzag a > a\, the 
chain becomes unstable with respect to the bending long- 
wave phonons, i.e., the bending instability of the chain 
takes place. Note that as regards the rest of (longitu- 
dinal and orientational) phonons, the chain keeps to be 
stable under the stretching a > ao). Thus, the bend- 
ing instability of zigzag-like chains does not admit the 
existence of the topological solitons of stretching. The 
maximum possible stretching of the chain can be defined 
as the ratio a±/l x - For the chain (• • -HF- • the maxi- 
mum stretching is 21%, while for the chain (• • HC1- • -)oo 
it is 29%. 

The similar bending instability of the chain is ob- 
served under stretching the trans-zigzag of the polyethy- 
lene macromolcculc ( — CH 2 — ) n . The analysis of the lin- 
ear dynamics of the planar zig zag of the chain within 
the model studied in Refs. [la, [l6[ has shown that the 
strained chain is stable under stretching (the longitudinal 
lattice spacing of the zigzag) a < ao, where the critical 
value ao = 1.745 A corresponds to the point of inflection 
of the effective potential of longitudinal stretching. How- 
ever, under stretching a > ao the chain becomes unstable 
with respect to the bending oscillations of the chain. Here 
the maximum possible stretching of the chain is 37% (in 
equilibrium the longitudinal lattice spacing of the zigzag 
is l x = 1.276 A). Here the bending instability is caused 
by the fact that the interaction of the second neighbors in 
the trans-zigzag occurs only because of the deformation 
of the valent angles CCC and it does not depend directly 
on the distance between them. The bending of the chain 
does not allow to break it without deforming the valent 
angles. 

Thus, in a zigzag-like polyethylene macromolecule, the 
energy of interaction of the second neighbors is of the 
same order as the energy of interaction of the first neigh- 
bors, but the formation of bistable ground states is im- 
possible due to the bending instability of the strongly 
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stretched chain. For the absence of this instability the 
angles of a polymer chain have to possess a sufficiently 
strong nonvalent interaction of the second neighbors. 
The molecular groups of polyethylene CH2 do not pos- 
sess the interaction of this type. However, the radi- 
cals CHR with sufficiently long chains can provide this 
interaction. The polyolefine macromolecules ( — CH2 — 
CHR — ) n : polypropylene (R=C3H7), polystyrene con- 
taining benzol rings in radicals (R=CgH5) and polyvinyl- 
carbazole have the required stability structure. In these 
macromolecules with strongly interacting side radical 
groups R under the strong stretching of the chain, the 
topological solitons can exist. 

The existence of the topological solitons of stretching 
can be expected also in the DNA double helix. Here 
the conformational interaction of neighboring sites of 
the sugar-phosphate lattice plays the role of the first- 
neighbor coupling, whereas the stacking interaction of 
the neighboring purine and pyrrolidine bases can be con- 
sidered as the second-neighbor interaction. The experi- 
ments on stretching a sins le DNA molecule @-(| exhibit 
the presence of a specific constant region (plateau) in the 
force-stretching diagram Q . At strength about 65 pN the 
non-typical behavior: the molecule becomes elongated at 
constant force up to 1.7 of its contour length. At further 
stretching the force again begins to grow. A similar be- 
havior under stretching exhibit also a-helices of protein 

Within our model this behavior will take place under a 
weak second- neighbor interaction, when K < K$. Here, 
in a certain interval of lengthening, the stretching occurs 
according to the two-phase scenario, when one part of 
the chain is found in a weakly and the other one in a 
strongly stretched state. In the force-stretching diagram 
a constant region (plateau) appears due to the stretching 
of the chain because of increasing only the portion of its 
strongly stretched part. Obviously, here the topological 
solitons describing the transition from the weakly to the 
strongly stretched phases of the chain have to exist. 

VI. CONCLUSIONS 

The study carried out in this paper shows under 
the stretching of molecular chains the conformational 



changes of these chains can occur that result in differ- 
ent ground states. The transition regions between these 
states can be described as topological solitons. In the 
simplest model with the nearest-neighbor interaction of 
the Morse-like type and the second-neighbor harmonic 
interaction, it is shown that under the chain stretching 
the ground state is realized as a regular configuration 
with alternating bonds ("long-short"). In this case, the 
chain can be found in two degenerate ground states ad- 
mitting the existence of topological solitons that describe 
the chain transition from the state "short-long bond" into 
the state "long-short bond". This situation is possible 
for the molecular chains with sufficiently strong interac- 
tion of the second neighbors. With weak interaction, the 
chain stretching leads to the appearance of one region 
with weakly and the other one with strongly stretched 
bonds. As a result of this non-uniform stretching, the 
presence of a broad plateau in the force-stretching dia- 
gram of DNA double helix and protein a-helix can be ex- 
plained. The boundary between the weakly and strongly 
stretched phases of the chain can also be described as a 
topological soliton. 

Finally, it should be noticed that the models with the 
second-neighbor coupling being responsible for the stabi- 
lization of homogeneous bistable ground states have been 
studied earlier [la [l7j - l20j . However, in these (diatomic) 
models, the effect of switching or controlling bistability 
by external forcing has not been considered. The bista- 
bility here has been attained intrinsically due to the re- 
pulsive interaction in the heavy-ion sublattice. 
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